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We demonstrate significant improvements of the spin coherence time of a dense ensemble of 
nitrogen-vacancy (NV) centers in diamond through optimized dynamical decoupling (DD). Cooling 
the sample down to 77 K suppresses longitudinal spin relaxation T\ effects and DD microwave 
pulses are used to increase the transverse coherence time T 2 from ~ 0.7 ms up to ~ 30 ms. We 
extend previous work of single-axis (CPMG) DD towards the preservation of arbitrary spin states. 

Following a theoretical and experimental characterization of pulse and detuning errors, we compare 
the performance of various DD protocols. We identify that the optimal control scheme for preserving 
an arbitrary spin state is a recursive protocol, the concatenated version of the XY8 pulse sequence. 

The improved spin coherence might have an immediate impact on improvements of the sensitivities 
of AC magnetometry. Moreover, the protocol can be used on denser diamond samples to increase 
coherence times up to NV-NV interaction time scales, a major step towards the creation of quantum 
collective NV spin states. 


In recent years, atomic defects in diamond have 
been the subject of a rapidly growing area of research. 
The most well-studied of these diamond defects is the 
nitrogen-vacancy (NV) color center, whose unique spin 
and optical properties make it a leading candidate plat¬ 
form for implementing magnetic sensors [H-fTcl as well 
as qubits, the building blocks for applications in quan¬ 
tum information. In particular, NV spin coherence times 
longer than a millisecond have been achieved in single 
NV centers at room temperature, either through careful 
engineering of a low spin impurity environment during 
diamond synthesis 11] or thr oug h application of pulsed 
and continuous [3, fl71 ] dynamical decoupling 
(DD) protocols. These long single NV spin coherence 
times have been instrumental in demonstrating very sen¬ 
sitive magnetic EM , electric jT|j , and thermal [li| mea¬ 
surements as well as high-fidelity quantum operations 

ES0. 

Achieving similarly long spin coherence times in en¬ 
sembles of NV centers can further improve magnetic sen¬ 
sitivity [1,@] and, moreover, may open up new avenues for 
studying many-body quantum entanglement. For exam¬ 
ple, achieving NV ensemble spin coherence times longer 
than the NV-NV interaction timescales within the ensem¬ 
ble could allow for the creation of non-classical spin states 


21-23]. Recently, NV ensemble spin coherence times 
up to ~ 600 ms have been demonstrated by perform¬ 
ing Carr-Purcell-Meiboom-Gill (CPMG) DD sequences 
at lower temperatures to reduce phonon-induced deco¬ 


herence [24]. The CPMG sequence preserves only a sin¬ 
gle spin component efficiently, however; experimentally, 
in the presence of pulse imperfections, the CPMG DD 
protocol cannot protect a general quantum state [2^2?} , 
as is necessary for applications in quantum information 
and sensing. To date, the preservation of arbitrary NV 
spin states has been considered only in a limited fashion, 
mostly at room temperatures and for single NV centers 


12-14J. However, no fundamental study yet considered 


the robustness of various DD protocols on NV ensembles. 
In this work, we perform a theoretical and experimental 
analysis of the performance of several DD protocols, in¬ 
cluding standard CPMG and XY-based pulse sequences 
as well as modifications thereon, and extract an optimized 
protocol for preserving a general NV ensemble state at 77 
K. We observe an extension of the arbitrary NV ensemble 
state from a coherence time ~ 0.7 ms of an Hahn-Echo 
measurement up to a coherence time ~ 30 ms, which 
is more than an order of magnitude improvement. Al¬ 
though higher coherence times were demostrated for pre¬ 
serving a specific spin state 24], in this work we funda¬ 
mentally study and optimize a DD protocol for preserv¬ 
ing an arbitrary state. 


The NV center is composed of a substitutional nitrogen 
atom (N) and a vacancy (V) on adjacent lattice sites in 
the diamond crystal. The electronic structure of the neg¬ 
atively charged NV center has a spin-triplet ground state, 
where the m s = ±1 sublevels experience a zero-field split¬ 
ting (~ 2.87 GHz) from the m s = 0 sublevel due to spin- 
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FIG. 1. (a) Energy levels of the negatively charged NV cen¬ 
ter, including the 14 N hyperfine splitting, A is the zero-field 
splitting, (b) Bloch sphere diagram illustrating the two main 
types of pulse imperfection: ej represents the deviation from 
an ideal rotation angle 7r, and n = ( n x ,n y ,n z ) is the actual 
rotation axis, which can deviate from k = ( k x , k y , 0). (c) Opti¬ 
cally detected magnetic resonance measurement of |0) | + 1) 
transition in an NV ensemble. Hyperfine interactions between 
the NV electronic and the 14 N nuclear spins form three NV 
resonances, and a strong static field ~ 300 G polarizes the 
14 N nuclear spins into the | — 1) spin state. 


spin interactions [Fig. [IJa)]. Application of an external 
static magnetic field along the NV symmetry axis Zee- 
man shifts the m s = ±1 levels and allows one to treat the 
m s = 0, +1 spin manifold (for example) as an effective 
two-level system. The NV spin state can be initialized 
in the m s = 0 state with off-resonant laser excitation, 
coherently manipulated with resonant microwave (MW) 
pulses, and read out optically via spin-state-dependent 
fluorescence intensity of the phonon sideband [lj. 

The NV spin bath environment is typically dominated 
by 13 C nuclear and N paramagnetic spin impurities, ran¬ 
domly distributed in the diamond crystal. These spin im¬ 
purities create different time-varying local magnetic fields 
at each NV spin, which can be approximated as a random 
local magnetic field that fluctuates on a timescale set by 
the mean interaction between spins in the bath. This 
random field induces dephasing of freely precessing NV 
spins on a timescale T£ |6, J, |28|, [29|. Dynamical decou¬ 
pling pulse sequences can suppress the effect of the spin 
bath noise and thus preserve the NV spin coherence up to 
a characteristic time 24, 29|. In the ideal case of per¬ 
fect pulses, various DD protocols (e.g., CPMG, XY, etc.) 
are equally effective at preserving an arbitrary NV en¬ 
semble spin state. Experimentally, however, off-resonant 
driving due to the NV hyperfine structure [301 ] and other 
pulse imperfections significantly affect the performance 
of individual DD protocols. In order to overcome these 
pulse imperfections, we optimize a DD protocol for an 
ensemble of NV spins. 

Figure 0(a) illustrates the general structure of the DD 


protocols explored in this work. In each protocol, (7r)- 
pulses about a rotation axis determined by the specific 
DD protocol are applied, with a free evolution interval 
of time 2r between them. In the regime where the pulse 
durations are short compared to the free evolution inter¬ 
val between adjacent pulses, each pulse can be expressed 
in terms of a spin rotation operator 2fij, 271 


U- k = exp {-Z7r(l + e k )[S ■ n]}. 


( 1 ) 


Equation m incorporates the two main types of pulse 
imperfection: e k represents the deviation from an ideal 
rotation angle 7 r, and h = ( n x ,n y ,n z ) is the actual ro¬ 
tation axis, which can deviate from k = (k x ,k y , 0) [Fig. 
Eb)]. Generally, imperfections in the rotation angle (e k ) 
may be caused by limitations in pulse timing resolution 
and amplitude stability of the MW field source, as well 
as static and MW field inhomogeneity over the measure¬ 
ment volume; and imperfections in the rotation axis may 
be caused by phase instability in the MW field source. 
In addition to general experimental pulse errors, the spe¬ 
cific physical system of the NV spin ensemble introduces 
additional pulse imperfections. Most notably, hyperfine 
interactions between the 14 N nuclear spin (I = 1) of the 
NV center and the NV electronic spin result in three 
transitions each separated by ~ 2.2 MHz in the, e.g., NV 
77is — 0 t-P +1 resonance [3l( [Fig. [l](c)]. 

The total evolution operator of a general DD sequence 
containing n (7r)-pulses can then be expressed as 


£/dd =U d {T)-U %n -U d {2r)-U hn i -U d (2T)-...-U d (2T)-U- ki -U d (T), 

( 2 ) 

where Ud is the free evolution operator. It is clear that 
without compensation for pulse imperfections in the spin 
rotation operators, accumulating errors will result in a se¬ 
vere loss of coherence even in the limit of free evolution 
time t —>• 0. First, we study the robustness of conven¬ 
tional CPMG and XY-based DD protocols, summarized 
in Figure [2] (b) (c), in order to determine which protocol 
is the most robust against pulse imperfections caused by 
general experimental limitations as well as those specific 
to NV ensembles. Realizing that enhanced robustness 
is necessary, we reduce the effects of the imperfections 
by optimizing experimental parameters (see detailed ex¬ 
perimental setup description below) and modify the ba¬ 
sic XY sequences by introducing pulses with additional 
phases [Fig. [2](d)] and concatenated cycles [Fig. [21(e)]. 
Similar DD protocol optimization has been performed in 
the past for phosphorus donors in silicon [261 ] and single 


NV centers mm mm- 


all ( 77 )- 
[Fig- Hb)]; 


In the conventional CPMG DD protocol 
pulses are applied along the same axis ( x) 
consequently, only coherence along one spin component 
is well-preserved. The XY family of DD protocols 0] 
applies pulses along two perpendicular axes ( x , y) in or¬ 
der to better preserve spin components along both axes 
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FIG. 2. Dynamical decoupling protocols. The directions of 
the arrows in the scheme represent the phases of the pulses. 
For each sequence, the free evolution time between pulses 2 t 
was swept to obtain a full coherence curve, (a) General DD 
scheme, (b) CPMG. (c) XY8. (d) KDD version of XY8: each 
(7r)-pulse from an XY8 sequence is replaced by five adjacent 
(7r)-pulses, with additional phases of (7r)6o° — (7r)oo — (7r)go° — 
(tt)o° — (tt)6o° , keeping a free evolution time of 2r between 
them, (e) Concatenated version of XY8: the first applied 
cycle (cycle 0) is a single conventional XY8. Each of the 
following cycles is constructed recursively from the previous 
ones: eight pulses of conventional XY8 are always applied, but 
between every two of them, the whole cycle from the previous 
iteration is applied. 


equally [Fig. [2(c)]. We also explored two DD proto¬ 
cols which introduce additional modifications on the ba¬ 
sic XY pulse sequences in order to improve its robustness 
against pulse errors. The first modification, the Knill 
Dynamical Decoupling (KDD) pulse sequence 12, tdj . 
introduces additional phases, thereby symmetrizing the 
XY-plane further and reducing the effects of pulse errors 
due to off-resonant driving and imperfect 7r-flips. In the 
KDD protocol, each (7r)-pulse in a basic XY sequence is 
replaced by five pulses with additional phases given by 
(tt) 60 ° - Mo° - M 90 ° - ( 7 r )o° - (tt) 60 °, where the 2r free 
evolution interval between adjacent (7r)-pulses timing is 
preserved [Fig. [2(d)]. The second modification employs 
concatenation, a recursive process in which every cycle 
is constructed from the previous cycles [Fig. [2(e)], and 
each level of concatenation corrects higher orders of pulse 
errors 


We performed measurements on an isotopically pure 
(99.99% 12 C) diamond sample with N concentration ~ 
2 x 10 17 cm -3 and NV concentration ~ 4 x 10 14 cm -3 
(Element Six), grown via chemical vapor deposition. The 
sample was placed in a continuous flow cryostat (Janis 
ST-500) and cooled with liquid nitrogen to 77 K, signifi¬ 


cantly reducing phonon-related decoherence to allow for 
NV spin coherence times 1 ms 24|, 37]. A 532-nrn laser 
optically excited an ensemble of ~ 10 4 NV centers within 
a ~ 25 /im 3 measurement volume, and the resulting flu¬ 
orescence was measured with a single photon counting 
module. A permanent magnet produced a static mag¬ 
netic field Bq ~ 300 G along the NV symmetry axis, 
Zeeman splitting the m s = ±1 spin sublevels. To coher¬ 
ently manipulate the NV ensemble spin state, we used a 
70-(im diameter wire to apply a MW field resonant with 
the m s = 0 <->• +1 transition. The spin rotation axes of 
the individual DD pulses were set through IQ modula¬ 
tion of the MW carrier signal from the signal generator 
(SRS SG384). 


As discussed previously, one of the sources of pulse im¬ 
perfections for NV centers is the hyperfine structure in 
the NV resonance spectrum; specifically, resonant driv¬ 
ing of one of the hyperfine transitions results in detuned 
driving of the other two, introducing both spin rotation 
angle and spin rotation axis errors. We mitigate these 
effects by: (i) applying a strong static magnetic field 
(~ 300 G) to polarize the 14 N nuclear spins [38( into one 
hyperfine state which we drive [Fig. [He)] and (ii) apply¬ 
ing a strong MW field to drive the NV transition with 
Rabi frequency (~ 15 MHz) much greater than the de¬ 
tuning due to NV hyperfine splitting (~ 2.2 MHz). Fur¬ 
thermore, we minimize general experimental pulse errors 
due to pulse timing and amplitude imperfections, MW 
carrier signal phase imperfections, and static and MW 
field inhomogeneities over the measurement volume |30}. 
We estimate that the pulse imperfections remaining af¬ 
ter this optimization are characterized by e^, ~ 0.15 and 
n z « 0.25. 


In order to determine how well each of the four DD 
protocols preserves a general NV ensemble spin state, we 
measure the NV spin coherence of two orthogonal initial 
spin components S x and S v . The S x spin component is 
prepared and measured by applying the initial and final 
(7r/2)-pulses about the y axis; likewise, the S y spin com¬ 
ponent is prepared and measured by applying the initial 
and final (7r/2)-pulses about the x axis. We first charac¬ 
terize the robustness of each DD protocol against pulse 
imperfections by measuring NV ensemble spin coherence 
in the short free evolution (i.e., decoherence-free) limit 
2 nr <C I 2 (while remaining in the regime of infinitely 
narrow MW pulses) and normalizing against the NV en¬ 
semble spin coherence of a 1-pulse Hahn-Echo measure¬ 
ment in the same limit. We plot the experimental re¬ 
sults in Figure [3(b) for each of the DD protocols as a 
function of number of pulses n, where a relative contrast 
of 1 corresponds to perfect preservation of NV ensemble 
spin coherence and relative contrast of 0 corresponds to a 
mixed state. Incorporating estimated pulse imperfection 
values into Equations m and ©, we also plot simulated 
relative contrast of each DD protocol as a function of 
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FIG. 3. Relative contrast in the decoherence-free limit (r -C 
of DD protocols as a function of number of pulses. For 
clarity purposes, the simulation is separated from the exper¬ 
imental results, (a) Simulation of the effect of non-ideal (re¬ 
pulses according to Equation ©. All XY8-based sequences 
performed similarly for initialization at S x and S y . (b) Exper¬ 
imental results. The relative contrast is determined via nor¬ 
malizing with a Hahn-Echo measurement in the decoherence- 
free limit. At the perpendicular axis, the contrast of XY- 
based sequences is similar, but the CPMG contrast vanishes 
completely, as demonstrated in the supplemental material [30| 


number of pulses [Fig. ©a)]. 

The CPMG protocol maintains the highest relative 
contrast for the spin component along the spin rota¬ 
tion axis of the DD pulses (S x ) but the lowest relative 
contrast for the spin component along the perpendicu¬ 
lar axis (S y ) 130] , as expected. The relative contrast of 
XY-based sequences is comparable for both spin compo¬ 
nents 30 but drops as the number of pulses increases, 
indicating that while the XY-based protocol is able to 
symmetrically compensate for pulse errors and thus pre¬ 
serve an arbitrary NV ensemble spin state, accumulating 
pulse errors due to imperfect compensation eventually 
limit the sequence to ~ 500 pulses. Within the XY fam¬ 
ily, we compared XY4, XY8, and XY16 pulse sequences 


34{ and found XY8 to show the best performance 


The KDD protocol, which introduces more spin rotation 
axes to further symmetrize pulse error compensation, and 
the concatenated protocol, which constructs the pulse se¬ 
quences recursively in order to correct for higher orders 
of pulse errors both improve upon the conventional XY8 


sequence, maintaining higher relative contrast for both 
spin components to > 500 pulses. Note that the mea¬ 
surements are in qualitative agreement with the simula¬ 
tions. Quantitavely, however, there is a disagreement, 
and the experimental results for the relative contrast are 
slightly lower than the simulation suggests. In particu¬ 
lar, the contrst of the concatenated XY8 protocol does 
not change with the number of pulses according to the 
simulation, which disagrees with the experimental data. 
This disagreement is likely caused by the interplay be¬ 
tween pulse errors and decoherence effects, which was 
not taken into account in the simulation and will be the 
subject of a future research. 
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FIG. 4. Experimental results of the coherence time of DD 
sequences as a function of the number of pulses, after initial¬ 
ization at S x - The results after initialization at S y are shown 
in the supplemental material [30 ] 


The measured NV ensemble spin coherence time is 
plotted as a function the number of pulses for each DD 
protocol in Figure [I] The CPMG, XY8, and concate¬ 
nated XY8 protocols all extend the NV spin coherence 
time as expected, given the nitrogen-impurity-dominated 
spin bath environment [2E|. However, the KDD proto¬ 
col is less effective at extending the NV spin coherence 
time; this underperformance is probably due to the fact 
that the phase difference between adjacent pulses in KDD 
(sometimes 60°) is smaller than in other sequences (90°), 
making phase errors more significant [30]. 

In conclusion, after optimizing experimental parame¬ 
ters to minimize pulse imperfections, we found the most 
robust DD protocol for preserving an arbitrary spin state 
in an NV ensemble system to be the concatenated XY8 
pulse sequence. By compensating for higher order pulse 
errors, the concatenated XY8 sequence maintains higher 
relative contrast than the conventional XY8 sequence and 
is expected to ultimately outperform the KDD sequence 
for larger numbers of pulses. Furthermore, the concate¬ 
nated XY8 sequence achieves longer NV ensemble spin 
coherence times than the KDD sequence. At 77 K, we 
measured an extension of the arbitrary spin state of an 
ensemble of ~ 10 4 NV centers by a factor of ~ 40 and up 
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to ~ 30 ms. 

The optimized DD protocol determined in this work 
may already have an immediate impact in improving the 
sensitivity of NV magnetometry [6| and, moreover, may 
be useful for quantum information applications. The 
sample in this work has nitrogen density ~ 2 x 10 1 ' 
cm -3 and NV density ~ 4 x 10 14 cm -3 , correspond¬ 
ing to N-to-NV conversion efficiency ~ 0.2% and typi¬ 
cal NV-NV interaction time ~ 150 ms. Using standard 
sample processing techniques, such as electron irradia¬ 
tion [j], to modestly improve the N-to-NV conversion ef¬ 
ficiency to ~ 1%, the concatenated XY8 pulse sequence 
can increase the NV ensemble spin coherence time to 
the NV-NV interaction time. In such a case, MREV- 
based techniques [39] can be applied to average out the 
NV-NV interactions and introduce effective Hamiltonians 
21-23], thereby creating self engineered quantum states 


(e.g. squeezed states) in NV ensemble systems. 
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